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Abstract 

The (planar) ERTBP describes the motion of a massless particle (a comet) under the 
gravitational field of two massive bodies (the primaries, say the Sun and Jupiter) revolving 
around their center of mass on elliptic orbits with some positive eccentricity. The aim of this 
paper is to show that there exist trajectories of motion such that their angular momentum 
performs arbitrary excursions in a large region. In particular, there exist diffusive trajectories, 
that is, with a large variation of angular momentum. 

The framework for proving this result consists on considering the motion close to the 
parabolic orbits of the Kepler problem between the comet and the Sun that takes place when 
the mass of Jupiter is zero. In other words, studying the so-called infinity manifold. Close 
to this manifold, it is possible to define a scattering map, which contains the map structure 
of the homoclinic trajectories to it. Since the inner dynamics inside the infinity manifold is 
trivial, two different scattering maps are used. The combination of these two scattering maps 
permits the design of the desired diffusive pseudo-orbits, which eventually give rise to true 
trajectories of the system with the help of shadowing techniques. 

Keywords: Elliptic Restricted Three Body problem, Arnold diffusion, splitting of separatrices, 
Melnikov integral. 


1 Main result and methodology 

The (planar) ERTBP describes the motion q of a massless particle (a comet) under the gravitational 
field of two massive bodies (the primaries, say the Sun and Jupiter) with mass ratio /i revolving 
around their center of mass on elliptic orbits with eccentricity e. In this paper we search for 
trajectories of motion which show a large variation of the angular momentum G = q x q. In 
other words, we search for global instability (“diffusion” is the term usually used) in the angular 
momentum of this problem. 

If the eccentricity vanishes, the primaries revolve along circular orbits, and such diffusion is not 
possible, since the (planar and circular) RTBP is governed by an autonomous Hamiltonian with 
two degrees of freedom. This is not the case for the ERTBP, which is a 2-I-I/2 degree-of-freedom 
Hamiltonian system with time-periodic Hamiltonian. Our main result is the following 

Theorem 1. There exist two constants C > 0, c > 0 and p* = p*(C,c) > 0 such that for any 
0 < e < c/C and 0 < p < p*, and for any two values of the angular momentum in the region 
C < Gl < G 2 < c/e, there exists a trajeetory of the ERTBP such that G(0) < Gi, G(T) > G 2 for 
some T > 0. 

*AD, AR, and TMS were partially supported by the Spanish MINECO-FEDER Grant MTM2012-31714 and the 
Catalan Grant 2014SGR504. 
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This result will be a consequence of Theorem where the large C and the small constant c 
are explicitly computed {C = 32, c = 1/8), and where it is also shown the existence of trajectories 
of motion such that their angular momentum performs arbitrary excursions along the region C < 
GI<G* 2 < cje. 

Let us recall related results about oscillatory motions and diffusion for the RTBP or the ERTBP. 
They hold close to a region when there is some kind of hyperbolicity in the Three Body Problem, 
like the Euler libration points |LMS85[ICZlllIDGR13] , collisions |Bol06) , the infinity |LS80IIXia931 
IXia92l IMoe071 IRob841 IMP941 IMS14| or near mean motion resonances [FGKR14] . Among these 
papers, two were very influential for our computations, namely [LS80) . where the Laplace’s method 
was used along special complex paths to compute several integrals, and [MP94| . which contains 
asymptotic formulas for a scattering map on the infinity manifold for large values of eG. Another 
one, |GMS12j . is very important for futute related work, since the proof of transversal manifolds 
of the infinity manifold is established for the RTBP for any /i S (0,1/2], 

Concerning the proof of our main result, let us first notice |LS801 [GMS12] that, for a non-zero 
mass parameter small enough and zero eccentricity, the RTBP is not integrable, although for large 
G its chaotic zones have a size which is exponentially small in G. This phenomenon adds a first 
difficulty in proving the global instability of the angular momentum G in the ERTBP for large 
values of G. 

The framework for proving our result consists on considering the motion close to the parabolic 
orbits of the Kepler problem that takes place when the mass parameter is zero. To this end we 
study the infinity manifold, which turns out to be an invariant object topologically equivalent to 
a normally hyperbolic invariant manifold (TNHIM). On this TNHIM, it is possible to define a 
scattering map, which contains the map structure of the homoclinic trajectories to the TNHIM. 
Unfortunately, the inner dynamics within the TNHIM is trivial, so it cannot be used combined with 
the scattering map to produce pseudo-orbits adequate for diffusion, and adds a second difficulty. 
Because of this, in this paper we introduce the use of two different scattering maps whose combi¬ 
nation produces the desired diffusive pseudo-orbits, which eventually give rise to true trajectories 
of the system with the help of the shadowing results given in |GMSI5| . 

The main issue to compute the two scattering maps consists on computing the Melnikov po¬ 
tential { 51) associated to the TNHIM. The main difficulty for its computation comes from the fact 
that its size is exponentially small in the momentum G, so it is necessary to perform very accurate 
estimates for its Fourier coefficients. Such computations are performed in Theorem and they 
involve a careful treatment of several Fourier expansions, as well as the computation of several 
integrals using Laplace’s method along adequate complex paths, playing both with the eccentric 
and the true anomaly. To guarantee the convergence of the Fourier series, we have to assume that 
G is large enough (G > G), and e small enough (Ge < c). Under these two assumptions, the 
dominant part of the Melnikov potential consists on four harmonics, from which it is possible to 
compute the existence of two functionally independent scattering maps. 

The combination of these two scattering maps permits the design of the desired diffusive pseudo¬ 
orbits, under the assumption of a mass parameter very small compared to eccentricity {0 < fi < p,*), 
see (801), which eventually give rise to true trajectories of the system with the help of shadowing 
techniques. 

It is worth noticing that since all the diffusive trajectories found in this paper shadow ellipses 
close to parabolas of the Kepler problem, that is, with a very large semi-major axis, their energy is 
close to zero, and the orientation of their semi-major axis only changes slightly at each revolution. 

The case of arbitrary eccentricity 0 < e < I and arbitrary parameter mass parameter 0 < /i < 1 
remains open in this paper. Indeed, the case eG « I involves the analysis of an infinite number of 
dominant Fourier coefficients of the Melnikov potential, whereas for the case eG > 1, the qualitative 
properties of the Melnikov function should be known without using its Fourier expansion. Larger 
values of the mass paratemer p than those considered in this paper involve improving the estimates 
of the error terms of the splitting of separatrices in complex domains, as is usual when the splitting 
of separatrices is exponentially small. The computation of the explicit trajectories from the pseudo¬ 
orbits found in this paper needs a suitable shadowing result given in [GMSIhj . which involves the 
translation to TNHIM of the usual shadowing techniques for NHIM. 

The plan of this paper is as follows. In Section we introduce the equations of the ERTBP, as 
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well as the McGehee coordinates to be used to study the motion close to infinity. In Section we 
recall the geometry of the Kepler problem, when the mass parameter vanishes, close to the infinity 
manifold and its associated separatrix. Next, in Section]^ we study the transversal intersection of 
the invariant manifolds for the ERTBP, as well as the scattering map associated, which depend on 
the Melnikov potential of the problem, whose concrete computation is deferred to Section The 
global instability is proven in Section]^ using the computation of the Melnikov potential, and is 
based on the computation of two different scattering maps, whose combination gives rise to the 
diffusive trajectories in the angular momentum. 


2 Setting of the problem 


If we fix a coordinate reference system with the origin at the center of mass and call qs and qj 
the position of the primaries, then under the classical assumptions regarding time units, distance 
and masses normalization, the motion g of a massless particle under Newton’s law of universal 
gravitation is given by 

Qs-q , 9j - 9 


X 9s - 9 , 9 j - 9 

Wi-Hi + ^- 

dt^ gs ^ 


9r I9J - 9l" 


where I — /i is the mass of the particle at qs and fi the mass of the particle at qj. Introducing the 
conjugate momentum p = dq/dt and the self-potential function 


Uffq,t;e) = 




1 - M _ 

\q-qs\ l9-9jr 


( 2 ) 


equation 0 can be rewritten as a 2-I-1/2 degree-of-freedom Hamiltonian system with time-periodic 
Hamiltonian 

(3) 


P 

Hffq,p,t-,e) = — - Uf,{q,t;e). 


In the (planar) ERTBP, the two primaries are assumed to be revolving around their center 
of mass on elliptic orbits with eccentricity e, unaffected by the motion q of the comet. In polar 
coordinates q = p(cosa,sina), the equations of motion of the primaries are 


qs= pr{cosf,sinf) gj =-(1 -/r)r(cos/, sin/). 


(4) 


By the first Kepler’s law the distance r between the primaries |Win4I[ p. 195] can be written as 
a function r = r(/, e) 

1 -e^ 

(5) 


T = 


1 -I- e cos / 

where / = f{t, e) is the so called true anomaly, which satisfies |Win4H p. 203] 

df {1 + ecosfff 


dt 


(l_e 2 ) 3/2 ■ 


( 6 ) 


Taking into account the expression Q for the motion of the primaries, we can write explicitly 
the denominators of the self-potential function ([^ 

l9-9sP = p'^-2p,rpcos{a-f) + (7) 
l9-9Jp = p^+ 2{1-p)rpcos{a-f) + {l-pfr'^. (8) 

We now perform a standard polar-canonical change of variables (g,p) '—> {p, a, Pp, Pa) 

f Pa Pa \ 

q = {pcosa, psina), p= (Ppcosa -^sina,PpSina- 5icQSQ,j 

to Hamiltonian 0. The equations of motion in the new coordinates are the associated to the 
Hamiltonian ^ 

H;{p,a,Pp,Pa,t-,e) = ^ ^ - U;{p,a,t-,e) (10) 
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with a self-potential U* 


From now on we will write 


U*{p, a, t] e) = Ufj_{pcos a, p sin a, t; e). 


G = F^, y = Pp 


so that Hamiltonian (10 1 becomes 




( 11 ) 


( 12 ) 


Remark 2. In the (planar) circular case e = 0 (RTBP), it is clear from equations ([^ and Q that 
r = 1 and f = t, and that the expressions for the distances 0 between the primaries depend on 
the time s and the angle a just through their difference a — t. As a consequence, U*{p, a, s; 0) 
as well as H*{p, a, y, G, s; 0) depend also on s and a just through the same difference a — t. This 
implies that the Jacobi constant H* -|- G is a first integral of system. 


2.1 McGehee coordinates 


To study the behavior of orbits near inhnity, we make the McGehee |McG73] non-canonical change 
of variables: 



(13) 


for X > 0, which brings the infinity p = oo to the origin a: = 0. In these McGehee coordinates, the 
equations associated to Hamiltonian (10) become 


dx 

dt 

1 3 

- 4 ^ y 

II 

lev- 

x^ dU^ 

4 dx 

(14a) 

da 


dG 

dU^ 


(14b) 

dt 

-x^G 

4 

dt 

da 



where the self-potential is given by 

Uf,{x,a,t-,e) = U*{2/x‘^,a,t;e) = %- (-—- + (15) 

2 \ as ajJ 


with 

l<7 - 9sP = 0's = 1 - cos(a - f) + ^p'^r'^x'^, 

l<? - 9 jP = o-j = 1 -f (1 - p)ra;^ cos(a - f) + ^{1 - pfr'^x'^. 


It is important to notice that the true anomaly / is present in these equations, so that the equation 
for / given in 0 should be added to have the complete description of the dynamics. 


2.1.1 Hamiltonian structure 


Under McGehee change of variables (131, the canonical form dp A dy + da A dG is transformed to 

4 

w =- irdx A dy + da A dG (16) 


which, on X > 0, is a (non-canonical) symplectic form. Therefore, expressing Hamiltonian (12) in 
McGehee coordinates 


y 

'Hf,{x,a,y,G,t-,e) = y 


x^G^ 

8 


Uf^{x,a,t;e), 


(17) 
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equations (14) can be written as 


dx 

x^ f 

dUA 

dy 

x3 f 

dt 

4 1 

. 9y ) 

dt 

4 1 

da 

dn^ 


dG 

577^ 

dt 

dG 


dt 

da 


57 ^ 

dx 


(18a) 

(18b) 


Equivalently, we can write equations (18) as dz/dt = {^,77^} in terms of the Poisson bracket 


’ ^ i \dx dy dy dx ) da dG dG da' 


(19) 


3 Geometry of the Kepler problem (/i = 0) 

3.1 The infinity manifold 



For fi = 0 and G > 0, Hamiltonian 0 becomes Duffing Hamilto¬ 


nian 


'Hoix,y,G) = 




X 

y 


( 20 ) 


and is a first integral, since the system is autonomous. Moreover, 
77o is also independent of e and a. Its associated equations are 


dx 


1 


dt 4 


^ _ b" 

dt 8 4 

dt 


(21a) 

(21b) 


where it is clear that G is a conserved quantity, which will be 
restricted to the case G > 0 from now on, that is, G € M+. 
The phase space, including the invariant locus a; = 0 is given by 


{x,a,y,G) G K>o x T x M x IR+. From equations (21) it is clear 
Figure 1: Level curves of 77o that 
in the (x > 0, y) plane, for 

fixed G > 0 £ao = {z = (x = 0, a, y, G) G K>o x T x K x M+j (22) 


is the set of equilibrium points of system ( |2l| ). Moreover, for any fixed a G T, G G M, 

A„,G = {(0,a,0,G)} 

is a parabolic equilibrium point, which is topologically equivalent to a saddle point, since it pos¬ 
sesses stable and unstable 1-dimensional invariant manifolds. The union of such points is the 
2-dimensional manifold of equilibrium points 

Acxd = Aq G- 

a,G 

As we will deal with a time-periodic Hamiltonian, it is natural to work in the extended phase 
space 

z = (z, s) = (x, a, y, G, s) G K>o x T x K x K_|_ x T 
just by writing s instead of t in the Hamiltonian and adding the equation 

ds ^ 
dt 
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to systems (181 and We write now the extended version of the invariant sets we have defined 
so far. For any a £ T, G G M, the set 


Aa,G = {5=(0,a,0,G,s),sGT} (23) 

is a 27r-periodic orbit with motion determined by ds/dt = 1. The union of such periodic orbits is 
the 3-dimensional invariant manifold (the infinity manifold) 

Aoo = [J Aq,,g = {(0, a, 0, G, s), (a, G, s) G T x K+ x T} ~ T x M+ x T, (24) 

a,G 

which is topologically equivalent to a normally hyperbolic invariant manifold (TNHIM). 
Parameterizing the points in Aqo by 

xo = xo(a, G, s) = (xo(a, G), s) = (0, a, 0, G, s) G Aoo — T x K+ x T 

the inner dynamics on Aoo is trivial, since it is given by the dynamics on each periodic orbit 

^t,o(xo) = (0, a, 0,G,s + t) = (xo(a, G), s -f t) = xo(a, G,s + t). (25) 


3.2 The scattering map 

In the region of the phase space with positive angular momentum G, let us now look at the 
homoclinic orbits to the previously introduced invariant objects. 

The equilibrium points Aq,g have stable and unstable 1-dimensional invariant manifolds 

7„.G = W^{K,g) = W\K,g) 

= \z = {x,a,y,G), T-Lo{x,y,G) =Q,a = a-G j -dx\, 

I J«o=o y J 

whereas the 2-dimensional manifold of equilibrium points Aoo has stable and unstable 3-dimensional 
invariant manifolds which coincide and are given by 

7 = VF"(Aoo) = IF"(Aoo) = {z = (x, a, y, G), noix, y, G) = 0}. 

The surface 

7a, G = W""(A„,g) = VF^(Aa,G) (26) 

= \ z = (x, a, y, G, s), s G T, T-Loix, y,G) = 0,a = a — G / —dx i (27) 

I J-Ho=o y J 

is a 2 -dimensional homoclinic manifold to the periodic orbit A^^g in extended phase space. 
The 4-dimensional stable and unstable manifolds of the infinity manifold Aoo coincide along the 
4-dimensional homoclinic invariant manifold (the separatrix), which is just the union of the homo¬ 
clinic surfaces ja.G- 


7 = VF"(Aoo) = IT'^(Aoo) = U 7a.G 

a.G 

= {z = (x, a, y, G, s), {a, G, s) G T x K+ x T, T-Loix, a, y, G) = 0} 


(28) 


Due to the presence of the factor —x^jd in front of equations (21), it is more convenient to 
parameterize the separatrix fa,G given in (26) by the solutions of the Hamiltonian flow contained 
in 'H.q = 0 in some time t satisfying (see |MP94j l 


dt _ 2G 
dr x"^ 


(29) 
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In this way, the homoclinic solution to the periodic orbit Aq, <3 of system (21 1 can be written as 


2 

(30a) 

ao{t] a, G) = a + TT + 2 arctanr 

(30b) 

yo{t,G) - 

(30c) 

Go(t;G) = G 

(30d) 

so{t; s) = s + t 

(30e) 


where a and G are free parameters and the relation between t and r is 




3 


)■ 


(31) 


which is equivalent to (291 on T-Lq. From the expressions above, we see that the convergence along 
the separatrix to the infinity manifold is power-like in r and t: 


^0: VOl 


a — ao — TT 


G 


2 


2 


T,t ^ ±00. 


(32) 


We now introduce the notation 


zo = Zo((T, a, G, s) = (zo(cr, a, G), s) 

= (a;o(o-; G), ao(o'; a, G), j/o(ct; G),G, s) G 7 (33) 

so that we can parameterize any surface 7 ct,G as 

7 a,G = {zo = zo(o-, a, G, s) = (zo(cr, a, G), s), cr G M, s G T}. 

and we can parameterize the 4-dimensional separatrix as 

7 = VF(Aoo) = {zq = zo(CT,a,G,s) = (zo(cr, a, G), s), cr G K,G G K+, (a,s) G T^}. (34) 

The motion on 7 is given by 

<^t.o(zo) = zo(cr ± t, a, G, s) = (zo((t ± t, a, G), s + t) (35) 

and by equations ( [^ , ( [M] ) the following asymptotic formula follows: 

0t,o(zo) - <)'t,o(xo) = (zo(cr ± t, a, G), s + t)- (xo(a, G), s + t) —0 (36) 

>-±oo 

The scattering map S describes the homoclinic orbits to the infinity manifold (dehned in 


(24)) to itself. Given x_,x+ G Aoo, we define 

>S'm(x-) x+ 


if there exists z* G W^(Aoo) H W^{A^) such that 

^ 0 fort->± 00 . (37) 

In the case fj, = 0 the asymptotic relation ( |36[ ) implies S'o(xo) = xq so that that the scattering map 
Sq : Aoo — > Aoo is the identity. 
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4 Invariant manifolds for the ERTBP (/i > 0) 

4.1 The infinity manifold 


In order to analyze the structure of system (18 1 , we will write T-L^ given in ( |17| ) as 

a, y, G, s; e) = Hoix, y, G) - yMAfj.{x, a, s; e) 


where we have written in (151 as 


(38) 


Uf_i{x, a, s; e) = Uo{x) + yAUf^{x, a, s; e) = — + y.AU^{x, a, s; e), 


(39) 


and we proceed to study the dynamics as a perturbation of the limit case y = 0. From (151, 


AUo{x,a, s;e) = Vim ALln{x,a,s]e) 
>-0 


^ + (^)Vcos(a-/)-^ 


[4 + x^r'^ + cos(a — /)] 


(40) 


where r = r{f,e) and / = f{s,e) are given, respectively, in 


For /r > 0, it is clear from equations (18 1 that the set £oo remains invariant and, therefore, so 
does the infinity manifold Aoo, being again a TNHIM, and all the periodic orbits Aq,_g also persist. 
The inner dynamics on Aoo is the same that in the case /r = 0, so that the parametrization Xq as 
well as its trivial dynamics (25) remain the same. 

4.2 The scattering map 

From [McG73| we know that W^(Aoo) and lF))(Aoo) exist for y small enough and are 4-dimensional 
in the extended phase space. The existence of a scattering map will depend on the transversal 
intersection between these two manifolds. 

Let us take an arbitrary Zq = (zo,s) = (zo(( 7 , a, G), s) S 7 as in (33). Now, we have to 
construct points in ^^(Aoo) and W))(Aoo) to measure the distance between them. It is clear from 
the definition of 7 that 

v= (VHo(zo),0) 

is orthogonal to 7 = VF“(Aoo) = IF®(Aoo) at Zq and then if the normal bundle to 7 is denoted by 

V(zo) = {zq -I- cr V, cr e K} 

we have that there exist unique points z^’'^ = ( 2 ;®’'^, s) such that 

{2r} = ^r(^oo) niv(zo). (41) 

The distance we want to compute between lF®(Aoo) and IF))(Aoo) is the signed magnitude given 
by 

d(io,M) =^ 0 ( 2 )))-'Ho(zp. (42) 

We now introduce the Melnikov potential (see [DGOOl iDLSOBj l 


/ OO 

AUo{xo{t] G), ao(t; a,G), s +1; e) dt 

-00 


(43) 


where AUq is defined in (40). Thanks to the asymptotic behavior (32) of the solutions along the 
separatrix and of the self potential close to the infinity manifold 

AUo{x, a, s; e) = 0{x'^) as x —>■ 00 

this integral is absolutely convergent, and will be computed in detail in Section 
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Proposition 3. Given (a, G, s) S T x K+ x T, assume that the function 

(7 S M I— C{a, G, s — a;e) G K (44) 

has a non-degenerate critical point a* = a-*{a, G, s; e). Then for 0 < fj, small enough, close to the 
point Zq = (zo(cr*, a, G), s) G 7 (see the parameterization in (33)), there exists a locally unique 
point 

z* =z*(a*,a,G,s;e,/r) G W^{A^) (h W;^{A^) 

of the form 

z* = zl + 0{^i). 

Also, there exist unique points x± = (0, a±, 0, G±, s) = (0, a, 0, G, s) -|- 0{fj.) G Aqo such that 


Moreover, we have 


't>tA^*)-(l)tA ^±)—>0 fort ^± 00 . 


dC 

G+ - G_ = G, s - a*{a, G, s; e)) + 0 (m"). 

oa 


(45) 

(46) 


Proof. From equation (33) we know that any point Zg G 7 have the form 

Zq = 2 o{a,a, G, s). 


As in (41), we consider 


= (zr,s)GtF«’'^(Ao.) niv(zo). 


and we are looking for Zg such that z® = z“ There must exist points x± = (z±,s) G Aoo such 


that 


^ t —>-±oo 


(47) 


moreover ~ 4't,o{^o) = 0{fi) for ±t > 0 (see [McG73) l. Since "Hg does not depend on 

time, by ( |38| ) and the chain rule we have that 

^^noAtAzA)) = {no,nAAtA^T)) = -AnoAuAAtA^Ay^e). 


Since "Hg = 0 in Ago, using (47) and the trivial dynamics on Aqo we obtain 

^±00 


p±oo 

■^0(2^’“) = -M / {'Hg,AG^}((^(,^(z^’");e)dt. 

Jo 


Taylor expanding in ^ and using the notation (35) 


/ OO 

{"Hg, AGg}(((it,g(zo); e) dt + 0{A) 

-00 

/OO 

{77g, AGg}(zg(tT -|-1, a, G), s -I-1; e) dt -I- 0(A)- 

-OO 


(48) 


On the other hand, from (43) 


and then 


so that 


/OO 

AGg(a;g(i^ — s; G), ag(iz — s; a, G), s; e) diz 

-OO 

dC f°° 

— (a, G, s;e) = - J {AGg,'Hg}(zg(iz - s,a,G),s;e) dv 
dC f°° 

— (a, G,s - cr;e) = / {"Hg, AGg}(zg(i/- s-|-cr, a, G), s; e) di/ 

os J_oo 

pOO 

= / {7tg,AGg}(zg(t,a,G),s-f t;e)dt 


(49) 
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and therefore, from (48) and (49) 


dC 

d{zo, ^l) = Ho(z“) - Ho(zP = G, s - cr; e) + 

For ^ small enough, it is clear by the implicit function theorem that a non degenerate critical value 
a* of the function (44) gives rise to a homoclinic point z* to Aoo where the manifolds W^(Aoo) 
and TF“(Aoo) intersect transversally and has the desired form z* = Zg + 0(/i). 

Consider now the solution of system (18) represented by By the fundamental theorem 

of calculus and (38) we have 


G+ - G_ = - 
= M 

= M 


dH, 

, da 
da 

dMAo 

da 




dAU^ 

da 




) dt 


{(j)t,o{zo);e)dt + 0{n'^) 

(zo(ct* +t,a,G),s + t] e) dt + 0{^j?) 


dC, 


= a*;e) + O(m^). 


□ 


Once we have found a critical point a* = a*(a,G, s]e) of (44) on a domain of (a,G,s), we can 
define the reduced Poincare function (see [DLS06| 'l 


C*{a, G; e) := G(a, G, s — u*; e) = G(a, G, s*; e) 


(50) 


with s* = s — a*. Note that the reduced Poincare function does not depend on the s chosen, since 
by Proposition!^ 

^ (£ (a,G, s - CT*(a, G, s;e);e)) = 0. 

Note also that if the function ( [44| ) in Proposition has different non degenerate critical points 
there will exist different scattering maps. 

The next proposition gives an approximation of the scattering map in the general case fi > 0 


Proposition 4. The associated scattering map (a+, G+, s+) = Sf^{a, G, s) for any non degenerate 

dC* 


critical point a* of the function defined in (44) is given by 

dC 


( OL, OL, \ 

a - G; e) + G + G; e) + sj 


where C* is the Poincare reduced function introduced in (50). 


Proof. By hypothesis we have a non degenerate critical point cr* of (44). By definition (50), 
Proposition gives 

G+ — G = fj,— —(a, G) + 0{fjf). 
da 

as well as G_ = G + O(^) to get the correspondence between G+ and G_ that were looking for. 


The companion equation to (46) 


dL* 


a+-a = G) + O(m^) 

is a direct consequence of the fact that the scattering map is symplectic. 

Indeed, this is a standard result for a scattering map associated to a NHIM, and is proven in 
|DLS081 Theorem 8]. For what concerns our scattering map defined on a TNHIM, the only differ¬ 
ence is that the stable contraction (expansion) along ^^’“(Aoo) is power-like (32) instead of expo¬ 
nential with respect to time. Therefore we only have to check that Proposition 10 in |DLS08j still 
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holds, namely that Area ((/)t_^(7?.)) —0 when t —0 for every 2-cell TZ in hh®(Aoo) parameterized 
by R : [0,1] x [0,1] —> fF^(Aoo) in such a way that R{ti,t 2 ) G ^^^(Aoo), -R(0,t2) G Aoo- But this 
is a direct consequence of the fact that the stable coordinates contract at least by C/^ (see (32)) 
and the coordinates along Aqo do not expand at all. □ 


Remark 5. In the (planar) circular case e = 0 (RTBP), a, s; e) depends on the time s and 

the angle a just through their difference a — s, see Remark]^ From 


da 


{x,a,s] 


0 ) 


dAU^ 

ds 


ix,a,s; 


0 ) 


one readily obtains 


and therefore 


dC dC 

— (a,G,s;0) = - — (a,G,s;0) 

dC dC 

— (a, G, s - cr*;e) = - —(a,G, s - cr*;0) = 0 


and consequently the reduced Poincare function C* does not depend on a, and G+ = G_ -|-0(^^). 

But indeed G+ = G_ in the circular case, since there exists the first integral provided by the 
Jacobi constant Gj = + G and as = 0 on A^o, G+ = G_. Therefore in the circular case 

there is no possibility to find diffusive orbits studying the intersection of W^{kao) and IF“(Aoo) 
since any scattering map preserves the angular momentum. 


5 Global diffusion in the ERTBP 


We have already the tools to derive the scattering maps to the infinity manifold Aqo, namely 
Proposition 1^ to find transversal homoclinic orbits to Ago and Proposition to give their expres¬ 
sions. Both of them rely on computations on the Melnikov potential C. Inserting in the Melnikov 


potential introduced in (431 the expression for AUq in (401 we get 


C{a,G,s]e) = f 

J —( 


-oo [ [4 -I- x^r"^ + 4a:Qr cos(ao — /)] 


1/2 


+ (y) rcos[ao - /) - y 


dt (51) 


where Xq and ap, coordinates of the homoclinic orbit defined in (301, are evaluated at t, whereas 
r and /, defined in (§ and ([^, are evaluated at s + t. 

To evaluate the above Melnikov potential, we will compute its Fourier coefficients with respect 
to the angular variables a, s. Since xq and r are even functions of t and / and ao are odd, C is an 
even function of the angular variables a, s: L{—a, G, —s; e) = L{a, G, s; e), and therefore L has a 
Fourier Cosine series with real coefficients Lq^k- 

£ = ^ ^ = Lo,o + 2 ^ To.fc cos fca -I- 2 ^ ^ Lq^k cos{qs + ka). (52) 

q^'Lk^’L k>l q>l k^Z. 


The concrete computation of the Fourier coefficients of the Melnikov potential (51) will be carried 
out in sectionFirst, some accurate bounds will be obtained: 
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Lemma 6. Let G > 32, q > 1, k > 2 and £ > 0. Then < Bq^ and |Lo,«| ^ where 

Bq^o = 


Bq^i = 2 V (1 + 

= 292V«eli-«lG-i/2e-«°'/3 

= 2 ® 2 '=e«(l + 


( 53 ) 


_fe = 252«+2'=e2«el'=-«lG'=-i/2e-9GG3 


■*q, — k 

i?n/ = 2822VG-2^-3 


Directly from this lemma, we first see that the harmonics Lq^^ are exponentially small in G for 
g > 1, so it will be convenient to split the Fourier expansion (521 as 


G = Go 


' • • — Go + Gi + G >2 


(54) 


where 


Go (a, G; e) = Lo,o + 2 ^ Go,fe cos fca, 

fc>i 

Gq(Q;, G, s; e) = 2^ cos(gs + fca), 9 > 1- 


(55) 


The function Go does not depend on the angle s and it contains the harmonics of G of order 0 in s, 
which are of finite order in terms of G, Gi the harmonics of first order, which are of order , 

and all the harmonics of Cq for q > 2 are much exponentially smaller in G than those of Gi, so we 
will estimate Gq and Gi and bound G> 2 - 


To this end, it will be necessary to sum the series in (55). From the bounds Bq k in (53) for 
the harmonics Lq^k we get the quotients 


Bq^k+i _ 2(1 + e) 


B, 


q,k 


G 2 


for k> 2, 


B, 


9 , — (A:+l) 
Bq, — k 


= 4eG for k> q, 


B, 


0,^+1 

Bq,e 


4e 

= ^ for e>o, 


(56) 


which indicate that, for fixed q, we will need at least the conditions G > \/2{l + e) and eG <1/4 
to ensure the convergence of the Fourier series. This is the main reason why we are going to restrict 
ourselves to the region G > G large enough and eG < c small enough along this paper to get the 
diffusive orbits. 


Among the harmonics Gg^fc of 0 order in s, by (56), the harmonic Gg^o appears to be the 
dominant one, but we will also estimate Gg i to get information about the variable a, and bound 


the rest of harmonics Gg^^ for k > 2. Among the harmonics of first order Gi^fc, again by (56), the 
five harmonics Li k for |fc| < 2 are the only candidates to be the dominant ones, but the quotients 


from (53) 


Bi.2 

Bi,-i 


(1 + e)^ 

8eG4 


^ 1,1 


(1 + e)^ 
8eG3 


l?i,o 

Bi,-i 


e 

G 


eG 

G 2 ’ 


(57) 


indicate that Gi__i and Gi _2 appear to be the two dominant harmonics of order 1. Summarizing, 
to compute the series (52) we estimate only the four harmonics Gg g, Gg i, Li _i and Gi _ 2 , and 
bound all the rest, providing the following result, whose proof will also be carried out in section]^ 


Theorem 7. For G > 32, eG < 1/8, the Melnikov potential (51) is given by 


G, sj e) — Gg(o, G] e) -t- Gi(o, G, s\ e) -t- G> 2 (o, G, sj e) 


(58) 


with 


Gg(Q;, G; e) — Gg g + Gg 1 cos 01 + Gg(Q!, Gj c) (59) 

Gi(a, G, s; e) = 2Gi,_i cos(s — a) + 2Li _2 cos(s — 2a) + Gi(a, G, s; e) (60) 
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where the four harmonics above are given by 


-^0,0 — + ^o,o) 

^ 157re_ ^ , 

■^04 - + 


Li^-i — 





Li,_ 2 = -3v^eG3/2e-GV3(i + 


( 61 ) 

(62) 


\atisfy 


-^0,0 

< 2^2G"^ + 22 49e' 

-^0,1 

< 2^^G-^ + 


< 221G"^ +2 49 6^ 

1-^1,-2I 

<2i^G-i + fe 

0 

l^ol < 2^4 e^G"^ 


|£l| < 2l9e-G^/3 \n-7/2 ^ g2g5/2 ^ g^-3/2- 

|/:>2| < 


(63) 

(64) 


The function Ci contains only harmonics of first order in s, so we can write it as a cosine 
function in s. Introducing 


p := =: 12eG2(l + Ep) 

Li^ — l 1 + £/ i ^_2 


(65) 


in the definition (|60|) of £i we can write 


Cl — 2Li _i j 


v,fcez 


Gi,fc 

Gi,-i 


= 2Li__i cos(s — a) — pcos(s — 2a) + 


= 2Li,_i5ft 1 - 


pe 


E 


^7^-1.-2 

L 


Ll,k 

Li-1 


cos{qs + ka) 


1;^ ^2(fc + l)a 


k^-l-2 


Li-i 


= 2LK_i3fJ = 2Li,_iBcos(s - a - 0) 


( 66 ) 


where B = B{a,G;e) > 0 and —6 = —0(a,G;e) S [—tt, tt) are the modulus and the argument of 
the complex expression 


1 - pe-“ + 




k^-1,-2 


^l.-l 


(67) 


Writing also in polar form the quotient of the sum in ([6^ by the parameter p introduced in (65) 


Ee-^^ := 

k^-1,-2 


Ll,k i(k+l)a _ 
pLp-i 


= - E 




1^1,k i(k+l)a 

Lr-2 


with E = E{a,G;e) > 0 and —(j) = —(j){a,G;e) S [—7r,7r), equation (67) for B and 9 reads now as 

Be-^^ = 1 - pe-*“ + pEe-^'^ (68) 
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or, equivalently, as the couple of real equations 


B cos 9 = 1— p cos a + pE cos (j) 
—B sin 9 = p sin a — pE sin cj). 


(69) 

(70) 


The function E = E{a, G; e) is small, since, by ( [M| , ( [6^ and ( [6^ , 

|fi| 


\E\< 


1 ^ 1 ,-2I \pLi^-i\ 


\£l\ ^(l + e)4G-7/2 ^ g2(^5/2 g(^-3/2 


1 


/tt 

37r(l + e)"’^\ 1 


3v^ 




(71) 


with an analogous bound for its derivative with respect to a. Writing equation (68) as 


Be-^^ = Be-^^ +pEe- 


one gets the explicit formulae for B and 9 


B = y/l — 2p cos a+p'^ = Y^(l — p)2 + 4psin^ a > 0, 

e (-7r,7r] 


(72) 


9 = —2 arctan 


psinof 


B + l—p cos a 


from which we see that B behaves like a distance to the point p = 1 and a = 0. The angle 9 is not 
well defined when _B = 0, but this happens only for a = 0 and p = 1, that is, for G ~ (12e)“^/^. 
A totally analogous property holds for B: 

Lemma 8. B{a, G;e) > 0 except for a = 0 and p = 1 + -2 Ei,k/Li,-i- 


Proof. For 5 = 0, equation (70) reads as 


sin a = f{a). 


(73) 


where /(a) = f{a,G;e) := Esincj) = E sin (j){a,G; e). Since + {df/daf^ < 1 due to ( [TI] ), 
there are exactly two simple solutions of equation (73) in the interval [—7r/2, 37r/2]; one is Og _(_ € 
(—7r/2,7r/2) obtained as a fixed point of the contraction a = arcsin (/(a, G; e)), and a second 
ttg _ G (7r/2, 37r/2) fixed point of the contraction a = n — arcsin (/(a, G; e)). Taking a closer look 
at equation (1^, we see that if a changes to —a, then —</), —9, B are solutions of (68) or, in other 


words, (j), 9 are odd functions of a and B even. Therefore a = 0, tt are the unique solutions of 


equation (70) for 5 = 0. Substituting a = 0,7r in (69) for 5 = 0, only a = 0 provides a positive 
p, which is then given by p = 1 +p5 = 1 + J2k^-i -2 ^i./c/Li,_i. □ 

We are now in position to find critical points of the function s 1 —G, s; e). To this end we 
will check that s 1 —>■ C{a, G, s; e) is indeed a cosine-like function, that is, with a non-degenerate 
maximum (minimum) and no other critical points. By Theorem The dominant part of the 
Melnikov potential C is given by Gg -I- Gi. By equation (58) and the bounds for the error term. 


for G large enough, the critical points in the variable s are well approximated by the critical 
points of the function Gg -I- Gi and therefore will be close tos — a — 0 = O,7r (mod 27r) thanks to 
expression (66). For this purpose, we introduce 


Cl = Gt(a,G;e) = 2Li,_i5 


(74) 


where 5 = B{a,G;e) is given in (67) and Gi^_i is the harmonic computed in (61). With this 
notation the function Gi can be written as a cosine function in s 


Gi(a, G, s; e) = Cl{a, G; e) cos(s — a — 9), 


(75) 
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and differentiating the Melnikov potential (58) with respecto to s we get 


dC 

ds 


—L\ sin(s — a — 9) + 


dC>2 

ds 


= 0 


sin(s — a — 9) 


1 dC>2 

C\ ds 


( 76 ) 


which is a equation of the form (73) for s — a — 9 instead of a and / = [dL> 2 /ds) jL\. Therefore, 
as long as B > {dC> 2 lds) /(2Li^_i), which by the estimate ( |^ for Li-i and the bound (64) for 

£>2 happens outside of a neighborhood of size O of the point 


(a = 0, G = G*) where G* « (12e) is such that p = 1 + ^ Li^k/Li-i, (77) 


there exist exactly two no-degenerate critical points Sg of the function s i—£(a, G, s; e). 

Let us recall now that the Melnikov function £, as well as its terms Cq are all expressed as 
Fourier Cosine series in the angles a and s, or equivalently, they are even functions of {a,s). 
Consequently, dCq/ds is an odd function of (a, s), and it is easy to check that each critical point 
s* is an odd function of a. Moreover, using the Fourier Sine expansion of dCq/ds, one sees that 
it s is a critical point of s i—>■ £(a, G, s; e), s -I- tt too, so sg,- = sg_+ -I- tt. We state all this in the 
following proposition. 


Proposition 9. Let £ be the Melnikov potential given in (58), G > 32 and eG < 1/8. Then, 
except for a neighborhood of size O of the point (a = 0,G = G*) given in (77), 

s I— C{a, G, s; e) is a cosine-like function, and its critical points are given by 


Sg^^ — Sg^_i_(Qf, G; e) — a -\- 9 -\- '®o,— — ^o.-i- Ttt — o, -\- 9 -\- tv -\- 

where 9 = 9{a, G; e) is given in (67) and = O 


From the proposition above we know that there exist Sg _ and Sg ^, non degenerate critical 
points of s !->■ £(a, G, s; e). Therefore, we can define two different reduced Poincare functions (50) 

£i(a,G;e)=£(a,G,sS,±;e) 

= £o(q!, G; e) ± C\{a, G; e) -I- E±{a, G; e). 


By the symmetry properties of £q, it turns out that each Cq{a, G; e) = £(a, G, tg ; e) is an even 
function of a. Moreover, since Sg _ = Sg _|_ -I- tt, one has that £* = (—1)'^£*, so we can write the 
reduced Poincare map as 


£i = £o ± £J + £; ± £* + £^ ± • • • (78) 

From the expression for the scattering map given in Proposition we can define two different 
scattering maps, namely 


a + £ QQ (o^j G] e) + G — fi (a, G; e) -I- sj . 


a£5 

da 


(79) 


These two scattering maps are different since they depend on the two reduced Poincare-Melnikov 
potentials C’^. From their expression (79), the scattering maps S± follow closely the level curves 
of the Hamiltonians More precisely, up to terms, S± is given by the time —/i map of 

the Hamiltonian flow of Hamiltonian The remainder will be negligible as long as 


ImI < 


dCl 


dCl 

dG 

1 

da 


Nevertheless, since we want to switch scattering maps, we will need to impose 

ImI < |£I| = 2|£i,_iH| = O (G-i/"e-«'/3) , 
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that is, ^ exponentially small with respect to G in the region C < G < c/e which a fortiori is 
satisfied for 


0 < n < fj.* =e 


(80) 


This is the relation between the eccentricity and the mass parameter that we need to guarantee 
that our main result holds. This kind of relation is typical in problems with exponential splitting, 
when the bound of the remainder, here is obtained through a direct application of the 

Melnikov method for the real system. To get better estimates for this remainder, one needs to 
bound this remainder for complex values of the parameter f or r of the parameterization (30) of 
the unperturbed separatrix. Such approach has recently been used for in the RTBP in |GMS12j 
and it is likely to work in the ERTBP, allowing us to consider any /i G (0,1/2], that is, imposing 
no restrictions on the mass parameter. 

We want to show now that the foliations of = constant are different, since this will imply 
that the scattering maps S± are different. Even more, we will design a mechanism in which we 
will determine the places in the plane (a, G) where we will change from one scattering map to 
the other, obtaining trajectories with increasing angular momentum G. To check that the level 
curves of C\ and C*_ are different, and indeed transversal, we only need to check that their Poisson 
bracket is zero. Since £!]_ and C*_ are even functions of a, their Poisson bracket {£!]_, £+} will be 
an odd function of a, so we already know that it will have a factor sin a. Using equation (78 1 we 
can write 


{£;,£;} = {£o + Ut + £; + •••, £o - Ut + u; - •••} = -2{£o, C\) + 

where £3 contains only Poisson brackets of odd order 

[9/2] 

odd>5 g—0 


Therefore, by the bounds ^ for the harmonics L,.*, the error term £3 = 0 is much 


and we now compute. 


exponentially smaller for large G than {Cq,C\}, which is O Te 

By splitting £ 0 , using (59), and £* = 2Li-iB, using (61), (72) and ( [7T| ), in their dominant 
and non-dominant parts 


£-0 — £-0 + £oi £1,-1 — + ^0,1); 


B = B + Eb, 


after a straightforward computation, we arrive at 




-2{£o,£l}-- 

where 


-2{£o,£*i} 

with 




d = 

^ 25 eG 5 

1-r COS a —— —— 


48 G 


—£\ STTpsina 
G4 

, 1 — cos a + p 

2G3 R2 


24eG 

G2 


and a small error term 


Ej = o(g-^ + eG-3 -b e^G^ -b pe^G'^ + p{eG + G-^))^G-^/^e-^'' 


5.1 Strategy for diffusion 

The previous lemma tells us that the level curves of £*^ and £*_ are transversal in the region 
G > 32 and eG <1/8, except for the three curves a = 0, a = it and d = 0 (which, by the way, are 
also transversal to any of these level curves, see figure]^. 
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e=0.06 



Figure 2: Level Sets of C\ {C*_) in Blue (Red) and d=0 in Green 


Thus, apart from these curves, at any point in the plane (a, G) the slopes dG/da of the level 
curves of and C*_ are different, and we are able to choose which level curve increases more the 
value of G, when both slopes are positive, or alternatively, to choose the level curve which decreases 
less the value of G, when both slopes are negative (see Figure]^. In the same way, we can find 
trajectories along which the angular momentum performs arbitrary excursions. More precisely, 
given an arbitrary finite sequence of values Gi, i = 1,... ,n we can find trajectories which satisfy 
G(T,)=G„i = l,...,n. 

Strictly speaking. This mechanism given by the application of scattering maps produce indeed 
pseudo-orbits, that is, heteroclinic connections between different periodic orbits in the infinity 
manifold which are commonly known as transition chains after Arnold’s pioneering work |Arn64) . 
The existence of true orbits of the system which follow closely these transition chains relies on 
shadowing methods, which are standard for partially hyperbolic periodic orbits (the so-called 
whiskered tori in the literature) lying on a normally hyperbolic invariant manifold (NHIM). Such 
shadowing methods are equally applicable in our case, where we have an infinity manifold Aoo which 
is only topologically equivalent to a NHIM (see |R.ob881 IR.ob841 IMoe021 IMoeOTl IGL061 IGLS14| 1 
and |GMS15) . 

With all these elements, we can finally state our main result 

Theorem 10. Let G\ < G 2 large enough and e > 0, /i > 0 small enough. More precisely 
32 < Gl < G 2 < l/(8e) and 0 < /r < /i* = . Then, for any finite sequence of values Gi G 

(Gi, G 2 ), i = 1, ■. ■ ,n, there exists a trajectory of the ERTBP such that G(Ti) = Gi, i = 1,... ,n 
for some 0 < Ti < In particular, for any two values Gi < G 2 G (Gj(,G 2 ), there exists a 

trajectory such that G(0) < Gi, and G(T) > G 2 for some time T > 0. 


17 






















diffusion along green lines 



Figure 3: Zone of diffusion: Level curves of £!j_ (£1) in blue (red) and diffusion trajectories in 
green. 


6 


Computation of the Melnikov potential: Proof of Theo¬ 
rem [7] 


The main difficulty to compute the Melnikov potential is that it is given by an integral (51) where 


the coordinates of the separatrix xq and oq are given implicitly (30) in terms of the time t through 
the variable t (311, whereas r and / are given in terms of s + t through the differential equation ([^ 
defining the true anomaly /. To evaluate the above Melnikov potential, we will compute its Fourier 


Cosine series (52) in the angles s,a. 


The next proposition gives formulas for its Fourier coefficients. To this we will consider the 
Fourier expansion of the functions: 




■»/(*) — 




n,m iqt 


(81) 


which can be found in |MP94) and [Wind 11 p. 204]. Using that r is an even function and that / is 
and odd function, one readily sees that the above coefficients are real and indeed they satisfy 

= c”’’" = (82) 

Once these coefficients c^’”* are introduced we can give explicit formulas for the Fourier coefficients 
of the Melnikov potential C. 

Proposition 11. The Melnikov potential given in © or in ( [5^ can be written as 

= with L, = ^L,,fce*'=“. (83) 

(jez feez 
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with 


/>1 

l>2 

L,,_i = 

l>2 

^9.fc = E'=9'"''’"''^(9,^-fc,0 for A: >2 

l>k 

Lg^-k = E 1,1-k) for k>2 


l>k 


and 


N{q, m, n) = 


2m+n /_l^2\ /—l/2\ f°° Qi‘}{'^+'^^/3) 


^ 2 m+ 2 n —1 \ ^ 


/_oo (r-'j)2™(T + z) 


2n 


dr 


Proof. We write Melnikov potential (511 as: 

C = Ci+ I rcos(ao-/) 


dt. 


(84a) 

(84b) 

(84c) 

(84d) 

(84e) 

(85) 

( 86 ) 


where 


can be written as 


A = 


' dt 


0 = / -^{l + ^rifit + s))e 


1/2 

-oo [4 + XqT'^ + AXqT cos(q;o — /)] 

(^1 + + s))e-*(““-^(*+"»^ 


*(“o-/(t+s)) 


- 1/2 


- 1/2 


Using 


we get 


where 


dt. 

(87) 


- 1/2 

I 


(i+z) 5 = E 

1=0 

-^1 = E E + E E 

/c>0 l^k k<.0 l<k 


^ (z -t) ( [^(/(^ + s))]"'-" dt; 0 < fc < Z 


ri — 

22 '- 


Sl = 


2-2l+k+l \ k — 




With these expressions is easy to see that Lq cancels out the last term in the integral (861 and 
that L\ + Szl cancels the cosine term, and so 


^ - E E + E *^-1 + E E + E E 


( 88 ) 


1>1 l >2 l <-2 

Now we perform the change of variable 

G3 


k>l l>k k<. — l l<k 


G'^ / T'"\ , G^ _ 2x , 

t=-(T + -), <ii=-(l+T )dr 
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introduced in pT|) , and we use the formulas for xq and Oq given in (30a I and (30b). In particular 

, 4 


we will use that: 


G2(l + r2) 


= 2GdT, 


e " = 


and the expansion in Fourier series given in (811, obtaining 




roo ^iq(T+T^/3) G^/2 




_ ika \ '' iq s 21 —k, — k 

- 


qI _ ikoL 

Oh. — e 


2-2i+fe ^_l/2V-l/2^ V- _2z+fc,_, r 


G-4'+2fe-l V -I \k- 


/ ' 


-oo {t ~ i) 2*(T + i)2(fc 0 


= e 


ika ^ ^ ^iq 2/+A;, —fc 


6 ^ C 


N{q, —I, k — 1) 


(89a) 

dr; I < k < —1 
(89b) 


substituting now equations (89a) and (89b I in the expansion ( 88 ) we get 


gGZ 1>1 

+ ^ ^ cE'’-'iV(< 7 , Z - 1,0 + E E -I, -I - 1 ) 

gGZ l>2 q^Z l<-2 

+ E E E l-kJ)+J2 E E -I, k-l) 


q^Z k>2 


l>k 


:/c <-2 


Kk 


(90) 


Now changing the indexes I —?> —I and k —>■ —k in the third and fifth terms one obtain the desired 
formulas for the Fourier coefficients Fg,/c- □ 


In view of proposition (11) and formulas (84), to compute the dominant part of the Melnikov 


potential and obtain effective bounds of the errors we will need to estimate the constants c( 


,n,m 

q 


defined in (81) and the integrals N{q^m,n) defined in (85) for q > 0 and only for indices m,n 
satisfying n > 0, m < n + 1. Alternatively to ([^, it will be very convenient to express the distance 
r between the primaries as 

r = 1 — ecosE (91) 


in terms of the eccentric anomaly E, given by the Kepler equation [Win411 p. 194] 

t = E — esmE. (92) 

This is done in the next three propositions. 

Proposition 12. Let n,m, q G Z, n,q > 0, m < n + 1. Then the Fourier coefficients defined 


in (81) satisfy 


I „ r29+"+ie9vd^el™-«l m > 0 

^ |(l + e)”''-i rn <—1 


''q I — 


Proof. In the integral formula for the Fourier coefficients 

c"’™ = — / r"e™^e-*«‘dt 

Stt In 


we change the variable of integration to the eccentric anomaly (92) {dt = rdE) to get 


n,m 


1 

= — (re*^)”'r"+4-™e-*«*dF;. 

27T Jo ^ ’ 


(93) 


(94) 
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To compute c"’"* from (94) we will use the identity (see [Win41l p. 202]) 


re 




— —e 
2a 


a = 


vTTe + vT^ 


which readily implies 


re*^ = a^e*'®' — e ■ 


,-i-E 


2 2 2 
2 ^ 2 ^ /;-^ 4 ^ 

“ +i? = '- “ = “ + 1 ^ 


4a^ 

^4 


= ^ae^E/2 _ ^^-^E/2Y 

2a 


= 1 - e^ 


a — 


16a4 


(95a) 

= ^/l^. (95b) 


To bound the integral (93) for m > 0 we will consider two different cases: 0 < q < m and 
0 < m < q. Let us first consider the case 0 < q < m. By the analyticity and periodicity of the 
integral we change the path of integration from ‘^{E) = 0 to ‘^E = ln( 2 a^/e): 

E = u + i\nf -M G [0, 27r] 


so that 


„iE _ iu-ln{2a^/e) _ ^ i 


and then, by (95a), (95b), and (92) 

re*-^ = ^e*“ — e + = e(cosu — 1 ) 


, e / e 2a^ 

r = 1 - - ^e™ H-e * 

2 \2a2 e 


^iu ^2^—iu 


= 1 - -r^e‘“‘ - a‘e 
4a2 


= 1 - 




e-*‘ = 


\4a2 
20 ^ 6 -*“ 


+„“) c«„+i (T, - «^) si„„ = 1 - 

exp \/l — e 2 cos u + i sin 


exp 


4a2 


e — a e 


2a^e 


therefore 


re^f\ = e(l — cosu) < 2 e 


jrj = \l{I — cosuf + (1 — e 2 ) sin^ u = — cos u) — sin^ u < 2 


2a^ 


\/l — e 2 cos 


2a^ 


= -exp ( — v 1 — cosu] < e 




Since 2a^ < 2, substituting these bounds in (94) we find directly the desired result for 0 < 9 < to. 
For the the case 0 < m < q we now perform the change of the integration variable through 

2 a 2 


E = V — i\u] - ) , r G [0, 27r] 


so that 


giB ^ g™+ln(2aVe) ^ 


and then, by ( |95a[ ), ( |95b[ ), and ( |92[ ) 

o3 


2n^ 

re*^ = -e™ — e - 


—e“™ = - 
e 


16a'^ 


cos V ——+ i \ a — 


16a^ 


sinu 


^cosr — ^(1 + cosu) + i'Jl — sinu^ 


r = 1 - —e 
4a2 


-™_aV’' = l-cosu 


sm V 


e-*‘ = 


ee 2 iv -iv 

^exp ae - ^e 


— i\/l — I 

exp ^\/1 — e 2 cos u + i sin uj 


ee 


2a2 
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therefore as before 


re*'^| < 2 , 


|r| < 2 


< —e^ 
- 202 ^^ 


Since 2a^ > 1, substituting these bounds in (941 we find the desired result for 0 < m < g. 

For TO < —1 we bound directly the integral (94) for E G [0, 27r]. Since je*-^! = |e“**| = 1 we have 

r27r 


K-l < ^ 


in+1 


dE 


by noticing that |r| < (1 + e) we conclude the proof of the bounds for the □ 

As we can see from equations (|84| the Fourier coefficients of the Melnikov potential C depend 


on the function N{q,m,n) defined in (85), so to bound (or to compute) these Fourier coefficients 
we need to bound (or to compute) N{q,m,n). 

Introducing 


I{q,m,n) = 


Loo (r-I)2’«(r + i)2 


-dr 


N{q, TO, n) can be written as 


N{q, TO, n) = 


Q2m+2n-l 


-l/2\ /-1/2 


/(g, TO, n). 


We will call 


Kt) = *(y + ^) 

the variable term in the exponencial of the integral, so that 


I{q,m,n) = J 


Qq^hir) 


-oo (t-z)2"*(t + j) 


2n 


dr. 


(96) 


(97) 


Since the integral /(g, to, n) involves an exponential, it will be useful the Laplace’s method |Erd56) 
of integration. In particular on a complex path with 3(/i(t)) = 0. So, let us define the path 


F = riUF2ur3ur4UF5 

where 0 < e < 1, t* is a point such that 9(/i(t*)) = 0 that will be hxed later, and 


(98) 


Fi = {r G C|3 (/i(t)) = 0} n {t G C|3?(t) < 3?(-rL} 

Fs = {r G C|3 (/i(t)) = 0} n {t G C|3?(t) > 3?( r*)} 

Fz = {r G C|3 (/i(t)) = 0} n {t G C|3?(-t*) < 3?(t) < 0} n {r G C||t - > ce} 

r 4 = {r G C|9(/i(r)) = 0} n {r G C|0 < sR(r) < SR( r*)} n {r G C||t - > ce} 

Fs = {r G C|9 (/i(t)) < 0} n {r G C||r — i| = ce}. (99) 



By means of the Cauchy-Goursat theorem and 
a limit argument, the integral /(g,TO,n), defined 
in (97) over the real axis, is equal to the one taken 


over the path F thinking of t as a complex number 
(see |LS80) '). In fact, by the same argument, its 
value does not depend on e. 

The positive branch of the hyperbola defined 
by S5(/i(r)) = 0 intersects the circumference of ra¬ 
dius e in two points that can be expressed as 
and —C'e defined by: 


c, = Fg n r 4 


- a = Fa n F, 


( 100 ) 


Figure 4: The complex path F 
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Since the integral over T does not depend on e, we will choose a particular value of e to bound 
I{q,m,n) and consequently N{q,m,n) defined in (85). Later on, in proposition 17 we will just 
compute the e-independent terms of this integral. 

It is not difficult to see that, if we define the function 


u{t) = h{i)-h{T) =-pr)= (r-if - ^(t - i)\ 


( 101 ) 


then 


1 + 


M(ri u r 2 ), u(r 4 u Fs) c 

Moreover, if r“ € Fi U r 2 then r+ = — f~ G F 4 U F 5 and 

u{t~) = 

On the other hand one can see that u is an increasing function while moving along Fi U F 2 or 
F 4 U F 5 in the direction of increasing imaginary part. Therefore u has two inverses in and 

r_. Before writing them down let us notice that the point defined in (100 1 can be written as 


Cr = i + ee 


iOp 


with G ( 0 , tt/2) 


and has the following expression in the coordinates u defined in ( 101 ) 


liCe) = \u{Ce) \ = e^\l- -ie*®'|= ke 


( 102 ) 


(103) 


with 


kg = |l-ie*^' 1 = 

I 3 I 


]J (^+ |sin6»£) +(^^cos0e'j 


> 1 , 


since by construction, 9^ G ( 0 , 7 r/ 2 ) and then 0 < sind^. 
Now, we can write the inverses of the function u 


:[u(C'e), -1-00) —F4 U F5 
u I— ^{u) + iri{u), 


T :[u(C'e),-1-00 )—S>FiUF2 
u I —> —C(u) + iri{u). 


(104) 


The change (101) is useful over Fi U F 2 and F 4 U F 5 , thus performing this change in (85), we 


have that for any £ > 0 
N{q,m,n) = 


Q2m+2n-l 




[-F)i,n(w)-^m.n(w)]e "dM-b(-f)e«''3 / /^,„(r)dr 




(105) 


where 


d — ? 2 

^m.n — 


m+n 


-l/2\ /-1/2 


m 


(u) = _ 

’ (t±(u) - z)2™+1(t±(w) -h f)2"+l 

(r) = 

J m,n \ ) 


^q^hir) 


(t — -I- f)2" ’ 


where h(r) is given in (96). Now several helpful lemmas follow. 


Lemma 13. Let m,n G Z, m,n > 0 and dm,n be defined by equation ( 106 ). Then 

Mm.nl < 0 - 1 / 22 ™+" if TO -I- n > 0 


(106) 

(107) 

(108) 


23 
















Proof. Let s G N, then 


- 1/2 

s 




5! 




1 


3 2s-1, 


1 

/ ix'* / 

1 , 

'' 1 \ 

< — 

(2- -) = I 

1 -) < lim 

1 -) 

- 2 « 

\ sJ V 

2 s/ s—> 00 ' 

^ 2 s/ 


2 s 

= p-i /2 


Next lemma gives information about the functions ni"^)- 


Lemma 14. The function „(it) defined in (107) has the expansion 


□ 


F±„(u) = (±v©)-^’"-^E©™’”(±V©)^' 

1=0 


where the coefhcients d™’" satisfy 


(109) 


d™’” = l/(2*)^"+\ |d7’”| < 1^- 

Consequently, the series (109) is convergent for \^/u\ < ©2/3. 
Proof. Let us introduce the function 


( 110 ) 


:= (±)a;^™+^F©„(a:© = ^d™’"(±a:)^ 

1=0 

which is well defined since m = is a good change of variables in M+ and has the two inverses 
X = ±a/u- To bound the coefficients d™’" we use Cauchy formula: 




Applying the change of variables 


|a;|=£ 


T^.nix) -1 

- —-—ax = -— 

27ri 








— 2m 


dx. 


= ±\l (r - i)^ - - 1)3 = ±(r - i)\l (1 - ©r - i)) = ±© 7 ©(a /2 - zr), ( 111 ) 


©3 


we obtain 

(±i)'d 7 ’’ 


_ AVlra-j 


= T 


(r- 1 ) 


2m.-j 


1 


3(1 — ir) 


27rl 


.,(0 — ? T ^ 2 

3^ ^ ^ (r-*)2™+i(r + z)2-+i2©3(2-zr)i 


fdr 


= T; 


1 *4^— 


/ 

J\r-iU 


dr 


2 27r3™ ^2 d|T-i|=p (t — l)t+3(T + l)2”(r + 21)^ 2 "* 

Now, taking p = 1 and using that |t + 1| > 1 and that 2 < |r + 21| < 4 we have 


K’”! < ( 3 


which is the desired bound. From this bound it is clear that the series defining „(a:) is convergent 
for I a: I < ©2/3 and therefore the one for „(u) is convergent for ©u < ©2/3. □ 
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From equation (1091 we have 


Fm,ni'^) = i±Vu) „(±v^), 

1=0 

where the regular part of the function „ (m) is given by 


5m,n(±V^) = (±V^) ^ d"'’'"iFVuY 

j—2m-\-l 


( 112 ) 


(113) 


and dj ’ are defined by equation ( |109[ ) and satisfy bounds ( |110 1 . 

Lemma 15. Let as in equation (113), 0 < /3 < 1 and 0 < y/u < /3-\/2/3. Then 

9 


|ffm,n(±VM)| < 


-)m—2 


1 -, 


Proof. It is clear from equation (113) that 


= ^<Q.2m+liFVuy- 


s =0 


Since by hypothesis 0 < y/u < /?-\/2/3 with /? < 1, we can apply lemma 14 to get 


l5m,«(±V^)l < q 


<,1 


s =0 


i: U (v^r 


EfI 


s=0 


2/3 = 


1-/3 


which proves the lemma. 

Next proposition gives a bound for N{q,m,n). 


□ 


Proposition 16. Let N{q,m,n) as defined in (105) for q > 0, m, n > 0, m + n > 0, C? > 1 . Then 

| 7 V( 9 ,TO,n)| < 2"+™+3e« 

Proof. We will bound the integrals of N{q, m, n) choosing 

e = C? > 1. 


We write down then, using (103) and that > 1, 


Fy.n{u)e du 


(Gd 

< |F± (w(C,))| 


< 


'G-^k^ 


l^m,n(w)|e du< / \F^Ju)\e ^ ^ du 


roc 3 Q'im+% 

/ e-«-“ du < - 3 - 

'g- 3 “ (2-(G-i))2’^+i gGS L 


e 2 


< 2G 


3m— I 


(114) 


It remains only the last integral of (|105|) where the integrand is given in (108) and the domain 

(115) 


Ta in (99). The path Fa can be parameterized by 

T = i + G“ 2 e‘^ with 9 G [ 6 >i, ^ 2 ] = [tt — 6^, 0 e] 


with 9^ given in ( 102 ). If we define 

h{9) = h{T{9)) = i 


r{9f 


+ , 
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a straightforward computation using (101) shows that 


h{0) = 

3 V 3i / 


3i 

and then, as G > 1: 

^g-fG3g-J(cos2S+iG-isin3e) ^ g-fG^ ^ J (1 +fG" § ) ^g-fG^gg^ 


Note that, by (115), over Fa we have that |t — z| = G 2<1 and therefore \t + i\ > 1, and we can 
bound the last integral of (105) using (116): 


^q^hir) 


Its (t - + i) 


2n 


dr 




gg^ii(6>) 


'G (r(0)-i)2™(T(6»)+z) 


2n 


de 


< 


|g9^ii(e)| 

4, |G-3/2|2- 


G"3 de < 


r»2 _3 


G-3m 


G"5 de < 7rG^™-3/2g-fG g,^ 


(117) 


From lemma 13 and the bounds (114) and (117), we can bound N{q,m,n) by equation (105) 
as follows 


\N{q,m,n)\ < + Tre‘^'^< 2™+’^+Ve-«^G 


Next proposition provides an asymptotic expression for N{q,m,n). 


•771 — 271— i 


□ 


Proposition 17. Let n + m > 0 and the constants d™’" be defined by equation (109) and d„_„ 
by equation (106). For q > 0, m,n > 0 and G > 1 we have 


N{q, m, n) = 


Q2771 + 271-1 




2ig' 


,s=0 


(2s- 1)!! 


im,n ^35—I j_ rpq . j^q 
^2m—2s^ ' ' -^rr 


m.n ' ^m.n 


where 


\T^J < 45 22^+2 . G-3 < 18q^-^G^"^-^. 

When s = 0 the factor l/(2s — 1)!! in the formula above should be replaced by 1. 

Proof. We proceed as in the proof of proposition [^changing the path of integration to the path F 
defined in (98) leading to the integral (|105[). The important fact is that the integral (105) does not 


depend on e. So, we will compute only the ^-independent terms of that integral. We will follow a 
series of lemmas leading to the proof of the statement. 

Lemma 18. Let 0 < e < 1 and u{Gg) be as in equation (103), defined by (107). For any 
e > 0 small enough we have, if G > 1: 

pCX> 2 2??1 pOO 2 

/ F±^„(M)e-«^“dw = V / e-«^“d™’”(±v^)-2™-i+Jdu + i? 

JuICt.) . 2_nJu(C,.) 


where the constants dj ’ are defined by equation ( |109[ ) and E satisfies 

\E\ < 45 22’"+2g-^ 

Proof. Let us take y/uf = f3y/2f3 with = —1 -|- ^^y^3 -I- -^11/2 ~ 0.79. A simple calculation 
using (101) shows that — i\ = \. By definition, for e > 0 small enough we have that 

0 < u(Ge) < u* < y/uf < y^2/3, so 


lu(Cc 


^m,n(w)e ^'^du = 


lu{C,) 


F^Ju)e-^^-du + E, 
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with 


which can be bounded as 


El = / F±,„(w)e-«^“dn, 


\Ei\ = 


E^Ju)e du 


< 


Q-q^u 


< 


2e-«— 


1 


|(r^('it) — z) 2 '"+i(t=*=(u) + z)^"+i 


1 


du 


qG^ |t=‘=(u^) — i|2m+i |r±(it*) + i|2™+i 
^ 22m+2^—3^—q^u^ ^ 22"i+2^—3 


(118) 


By lemma 14 and equation (112) we have 

/.M* „ 2m 


= V / d™’”e-«^“(±yu)-2—1+^dn + E 2 

J-(Ce) 


where 


i;2 = 


lu{C,) 


5mm(±V^)e = “dw- 


Using that = I3^j2j^, by lemma 15 we have that, for any £ > 0 small enough, 

^ /■“* , g3 roo 

\E2\< |g±„(±v^)|e-«^“dri<9— / e-«^“dn 

Ju{Ce) ^ — P Jo 

om—1 om—1 

< 9 _ G~^ < 9 - _ G~^ 


Finally, 


4Ce) 


d7’”e-«^“(±v^)-2™-i+^dM= / d7’"e-«^“(±v^)-2™-i+Jdu + i33(j) (119) 


dCe 


where 


^00 3 

E:i{j) = - / d7’"e-«'^“(±v/^)-2™-i+^dn. 

J U* 

Let us bound E^{j) using the inequalities of Lemma 


14 


|^^3(j)l < Mr’”|(v^)-'™-^+' / < |d™’"|(V^)-2™-i+^'2e-«'^“ 


i G 


-3 


<2|d7’”|(/?V2/3; 


— 2m—1+j 


G-^ <2{ ^ 


-2m-l+j 


G 


-3 


_ g 2 ”i—lyj—2m—l+j ^—3 

then, calling E 3 = J 2 ^=i Esij), we have 


2m 


IU 3 I < 9 2”*-iG-3^/3-2™-i+i < 9 2 '"-ig-3^ 


j=o 


-P 


Now the lemma is proven using that Xjfi < \f2 and 

\E\ = \Ei + E 2 + E^l 


□ 


27 















The next lemma is a straightforward application of the last one. 


Lemma 19. Let 0 < e < 1 and u{Cs) be as in equation (103), 
small enough we have, if G > 1: 


defined by (1071 for any e > 0 



[^m,n(u) - F^^^iu)]e ^du = 


^ poo 3 

s^oJu(C,) 


where E is the same as in lemma 18 


Proof. By lemma 18 we have 


/‘OO 2 2171 ^ 

■Jn{Ce) ^O-Jy-iCe) 

then the non trivial terms in the sum are given when —2m — 1 + j = —2s — 1 with s = 0,..., m. 
This observation proves the lemma. □ 


Lemma 20. Let 0 < £ < 1 and u{Cg) be as in equation (105). Then the e-independent term of 


lu{Ce) 


e « 2 “dS;- 2 s(Vw) irfu 


IS 


(s-b 1)! 


(-l)^ 2 «+^( 2 s + r(l/ 2 ) 


(2s-b 2)! 

Proof. By equation (103) we know that u{Ce) = 0{e^) and then the following definitions make 
sense, calling S = 


poo 

4,(e)=/ 

Ju{Ce) 

fp,s{s) = «(G,)P-(2«+i)/2e-9-5-(Ce), 
Using this notation and integrating by parts we have 




q6 


r°° '\P-(2s-l-l)/2 p-(j<5«(Ce) 

Ju(Cs) p- s-l/2 

iqdlp^sie) - fpA^)) 


p — s — 1/2 
1 


and also 


p — s — 1/2 

pOO 3 

/ = dAA-2shAe)- 

Ju(Ce) 


Now, in the case where s > 0, using equation (120) s times we get 

{qsy 


(-S - 5 -b l)(-s - i -b 2) • • • (-i) 
The e-independent term of Io,si£) is given by 

_ (AY _,. r / ^ _ 

(-S - i -b l)(-s - i -b 2) • • • (-i) e™ 


YAY - E 


( 120 ) 

( 121 ) 


{qSY Vp,s(e) 

^ (-S- i -bl)---(-s- i -bp)' 


(qsy 


(—s — I + 1)(—s ~ 5 + 2) ■ • • (—|) AA 

(v^)2-l 


r(i/ 2 ) 


(_s-i + l)(-s-i+2)...(-i) 


r(i/ 2 ). 
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Then the e-independent term of the integral in equation (121) is 


(-S 


i+l)(_,_i+2)...(-l) 


r(i/ 2 ) 


when s > 0. 

In the same way, we have that the e-independent term of 


/f 


/■OO 3 

,o(e)=/ 

Ju{Ce) 


is ^ r(l/2). Therefore the lemma is proved if we notice that 

( — s—--|-l)( — s—--|-2)---( — -) = -^^(25 — l)(2s — 3) • • • (1) 

(-1)" (2s-hl)!! (-1)" (2s-f 2)! 

“ 2« 2s-fl “ 22^+i(2s-f 1) (s-hl)! 

where we have used that 


(2s-f 1)!! = 


(2s-h 2)! 


2'<(s-hl)! ■ 

This expression allow us to write the cases s > 0 and s = 0 in one equation which completes the 
proof. 


□ 


Next Lemma is a straightforward application of lemmas [T^ and 


Lemma 21. Let u(C'e) given in equation (103) and defined by (107), then the e-independent 
terms of 

poo 3 

/ [-Pm.n(w) - F- „(u)]e-«^“dM 

Ju{Ce) 


are given by 


m / _|_ 1 \| 

^(_l)^2^+i(2s + r(l/2) + 2E 

B=0 


where E is the same as in lemma 18 


Lemma 22. Let „ be defined in equation (108), then 

1 f-qG^ 


1 / 

Res(/^_„(r),i) = 2ie-«‘^"/3 - f 

1^0 ^ 


^2m-l-2Z 


( 122 ) 


Proof. We use the definition of ^ given in (108), with /i(r) given in (96). Now, using (101), 

/i(t) = —2/3 — (r — -I- i(r — i)^/3 

and we have, taking any J > 0 small enough, 

1 


Res(/^_„(T),i) = 7^. f f^,ni'r)dT = 7^. [ 
2tti J\r-i\=s 27rz Ji, 


Q<l^h{T) 


27rz Jir-i\=6 {t - + i) 


2n 


dr. 


We use again one of the changes (111), for instance 


X = \/h{i) - h{T) = J W‘ 2 - - it), 

v3 
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obtaining 


R-es (/^,„(r),i) 


g-qG^S r ^-qG^x‘^/2 

TT J\x\=s {t+{x) - i)2""+i(r+(a::) + i)2n+i 

2ie-««'/3Res (xF”’'"(a;2)e-«^o) . 


Now we can use the Taylor expansion of the function ^ 

expansion of e“'?G^a :72 _ J^i^q ^ {—qG^x^/2')'‘ to obtain the desired formula. 

□ 


From this lemma one and the bounds for d™’" given in (110), one has 
|Res (/^,„(r),0| < 3 ^ 


1=0 


< 3 2’”+ie"'5"^"/^ 


om+l „m—l/^3m—3 
^ 9 <-T .-oG=‘, 


</G73 


(123) 


Now we can prove proposition EH N{q,m,n) is given in ( |105| ), and since it does not depend 
on £ we can apply lemmas 21 and 22 and bound (123) to obtain 


N{q, m, n) = 


d 

^m,n~ 

Q2m-\-2n-l 


m z' I 1 V 

^ (_l)«2«+i (2s + 1) ^ ® r(l/2) + 


.s=0 


where 


and by lemma 


Rrn.n = (“^e® " / fLni'^)dr 


By lemma 22 


\R 


q I < 

m.nl — 


\T^ J = 2F; < 45 22'"+2 . G-^ 

2^+1^771—1 


3 m—2 


-G 


3m—3 


m—l/^3m—3 


< 18g™-'G' 


Using that 2®+^(s + l)!(2s + 1)!! = (2s + 2)! to show that 

(2s + l)(s + l)! 1 


(2s+ 2)! 


2"+i(2s- I)!!’ 


completes the proof of the proposition Due to the fact that the right hand side of this last 
expression is not defined when s = 0 but the left hand side is and is equal to one, we need to point 
out that when s = 0, the term l/(2s — 1)!! in the final formula should be replaced by 1. □ 

The proof of Theorem [7] will be done constructively through the following series of lemmas and 
propositions. 

Let us first compute some coefficients c”’"*, more precisely c+, c^’^, Cg’° and Cq’^ 


Lemma 23. Let 


be defined by (81). Then 


3,1 


— 1 + Qi 


c+ = -3e + Q 2 , 


2,0 


— 1 + Q 3 ; 


3.1 _ J I ^ 

Cg ~ ~ 2 ^ Qi, 


with 


\Q^\<98e^, z=l,2,3,4. 
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Proof. From its definition given in (811 and using the change of variable t = E — e sin if we have 




2.0 1 f 3,lf 




3._ 

Cn = 


From equations (95) we have 


cjd = ^ - e + fwe-*^](l - ecosF;)^e-*^e*"""^d£; 


27r 


27r 


Aa? 


r2Tr 


[aV^ - e + (1 - ecosE)e-^^e^^^'^^dE 


2,0 


1 

= — / (1 — e COS if )^c?if 

Stt Jq 


3,1 
Cn = 


4 


[a^e*^ — e H-^e *^](1 — ecosif)^(iif. 

4a2 ^' 


In what follows we will use (see (95b)) that 

1 


0 < e < 1, - < < 1, + -^ = 1. 

2 4a2 


To bound we use equation (124). It is easy to see that 

a^e*^ — e + -^e“*'® = e*^ — e + ifi, 

4a^ 

(I — ecosif)^ = 1 — Secosif + E 2 , = I + ie sin if + ifa 


where 


(124) 

(125) 

(126) 

(127) 

(128) 


(129a) 

(129b) 


El = (a^ - l)e*-® + 


e 

4a^ 


1 (ie sin if ) 1 +^ 


satisfy 


e e 


E 2 = Se^ cos^ F; - cos^ E, £3 = - V 2- od 

0 + 2 )! 


^^l|<^ + ^=e^ |£;2|<4e^ < e^-< 2 e^ 


Using equations (129), we have from equation (124) that is the Fourier coefficient of order 1 
of the function 

(e®'® — e + ifi)(l — Secosif + if 2 )(l + iesinif + E^) = 

e*^ — e — 3e cos ife*^ + ie sin ife*^ + Qi (if) 


where 


Qi{E) =Ei — 3e^ cos if — 3eifi cos if + if 2 (e*^ — e + ifi) — sin if — 3ie^ cos if sin ife*'® 
— 3ie^ cos E sin E — 3ie^ sin E cos EE 2 + ie sin EE 2 (e*^ — e + if 1 ) 

+ E^ie^^ — eP El — Secosife*^ — 3e^ cos if — 3eifi cos if + if 2 (e*^ — e + ifi)), 


which implies that, up to order one in e, the Fourier coefficient is exactly 1. From the bounds 
for ifi, E 2 and E 3 we find |Qi(if)| < 98e^, which implies the lemma for cf'^. 

From equation (129), it is easy to see that 


a^e*-® - e 




4a2 


= [e*^ - 


.£;il^=e 2 *^- 2 ee*-® 


E 4 
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where 


Ei = + 2 Ei{e^^ - e) + Ef 

can be bounded, in regard of equations p 8 t and the bound for Ei, as 

{Eil <e^ + 2e^(l + e) + < 6e^. 


Using equation (129), we see from equation (125) that c-^ 
the function 


is the Fourier coefficient of order 1 of 


(e 2 i£ _ 2ee*^ + Ei){\ — ecos£')(l + iesinE + E^) = 

_ g Ee^"^^ — 2ee*^ + ie sin Ue^*^ + Q 2 (-£') 


where 


Q 2 {E) = 2e^ cos Ee^^ + E 4 — eE^ cos E 

= ie smE{—e cos Ee^^^ — 2ee*^ + 2e^ cos Ue*^ + E 4 — eE^ cos E) 
= i? 3 (e^*^ — e cos Ee^^^ — 2ee*^ + 2e^ cos Ee^^ + E 4 — eE 4 cos E). 


From the above expressions we conclude that, up to order one in e, the Fourier coefficient is 
exactly —3e, and from the bounds for E 4 and E^^ we find that \Q 2 {E)\ < 50e^ which implies the 


lemma for Ci 


We compute Cq’° using equation (126), as well as equation ( 129b[) to get 


^27r 


Cq’ = — / {I - 2iecosE + E2)dE = I E Qs 
27r ./n 


with 


o2ir 


Q, = - E2dE 


and we have immediately, using the bound for E 2 , that IQ 3 I < 4e^, the desired result for Cq . 
Finally, we compute Cq’^ using equation (127), as well as equations (129b) 


1 


= — / (e*-^-e + £;l)(l-3ecos^; + F;2)d^^■ 

27^ Jo 

Now, we want to find, up to order e the Fourier coefficient of order zero of the function 
(e*^ — e + Ui)(l — 3ecosU + i? 2 ) = e*^ — 3ee*^ cosE — e + E^, 


where 

E 5 = £’ 26 *^ + 3e^ cos E — eE 2 + Ei — 3eEi cos E + E 2 E 1 , 


from where we find 


3,1 



Qi 


with 


and using the bounds for £2 and £ 1 , 


1 - 

Q4 = — y E^dE, 
we find IQ 4 I < 19e^. 


□ 
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Lemma 24. Let G < 32. It q,k G N, q > 1, k > 2, then the Fourier coefficients of the Melnikov 


potential (52) verify the following bounds: 


|L,, o| <29(2e2)«e«G-3/2e-«GV3 

\Lq, i| < 2^e«(l + e)^ 

< 2^ (2e2)^ 

\Lq, k\ < 2®2'=e«(l + ef 

|L,,_fe| < 2^2^'= {2ey 

Proof. From equations (84) by propositions and we have 

\Lq, o| < 24e«e-«®'/3(2ee^^i^)‘?G-i/2 ^(2^0-1)' 


/>1 


\Lq, i| < 2 ^e‘^e-’>° ^((1 + ef2'^G-^y 

l>2 

\Lq,-i\ < e'?e-«‘^'/3e«^"i^2«el^-«lG3/2^(2^G-i)' 

l>2 

\Lq, k\ < 242-*^e«e-«‘=''/3(l + e)-'=G-*^-^/2 e)^2^G-y 

l>k 

\Lq_k\ < 242-2'=e«e-«‘^"/3e«'^^^2«el'=-«lG2'=-i/2^(24G-i)' 


l>k 


since by hypothesis 2^/G <1/2. Since all these series converge we have proven the lemma using 
that 0 < e < 1. □ 

Lemma 25. If q € N, q > 2. Assume G > 32, eG < 1/8, Then for q > 2 

\Lq\ < ^ 2i3e-«'^'/3(e223G)«G-i/^ (131) 


Proof. From lemma (24) we have: 


\Lq,k\ < \Lq,o\ + \Lq,l\ + iLq.-ll + ^(|ig,/c| + \Lq-k\) 


feGZ 


fc >2 


< e 


-qG^/3 


2®2'?e«e2«G-3/2 + 2^(1 + e)^e‘^G-’^/^ + 2^2‘^e‘>-'^e^‘>G-^/'^ 


k>2 

3/ 


< e-«^ 


2i°2«e«-ie29G-i/2 + 2^(1 + e)VG-'^/2 + 2^e‘^G-^/^ 


k ^2 


q-1 


+ 2^%^‘iG-^lH^e.^ ^(4Ge-i)'= + 2^%^‘^G-^l'^e-n^ Y^^'teGf 

k—q 


fc=2 
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Using that eG < 1/8 




2i029g9-ie29c'-i/2 




2-25g5-iG-9 2-3‘?e"2'JG"^"« 

< 2i3e-«‘^'/3(e223G)«G-i/2 


+ 2-39e-29G-^-« + G-^e + 1 

which concludes the proof. □ 

The Melnikov potential £ (51) has a Fourier Cosine series (52) which can be split with respect 
to the variable s as £ = £o +£i +£2 + • • •, like in ( 54p5 ), as well as a complex Fourier series (83) 
£ = Both series are related through Cq = Lq and £, = 2^ \e^‘i‘‘Lq] for g > 1. In 

the next lemma we see rather easily that the terms £>2 = £2 +£3 + £4 + • • • of second order with 
respect to s satisfy a very exponentially small bound in G. 

Lemma 26. For G > 32, eG < 1/8, £> 2 (a, G, s; e) = 2 EE £q,fc cos(gs + ka) is bounded 

^>2 fc^Z 


as 


|£> 2 (a,G,s;e)| < . 


Proof. By lemma [25| 


|£>2|<233g-i/2^ 

q>2 


-GV3^2^3g 


1 <? 


< 220e4G3/2e-2GV3 


where the last bound holds if 

g-2G3/3g223^ < 1/2 

which is true for every G > 32. Now, using that e < 4 we get the result. 

The next step provides an asymptotic formula for £1 = 23? {e®®£i}. 

Lemma 27. For G > 32 and eG < 1/8 we have the following formula for Li (83) 

3? 


(132) 

□ 


where 


{e*Ui} = 3? {e** ((c?’^iV(l, 2,1) + £ 3 ) 6"*“ + (c^’'A^(1, 2,0) + £i_.2) e'^^ + £ 1 ) } (133) 

|£i(a,G;e)| < 233e-^'/3 ^q-3/2 ^ q-7/2 ^ ^ 2 ^ 5/2 

|£ 3 (a,G;e)| < 22V^'/3g-3/2 
|£ 4 (a,G;e)| < . 


Proof. From equation (83), we have that 

£1 = Li.o + ^(£i,fce*'=“ + 
fc>l 

= £i._ie-“ + £i,_2e-2“ + ^ £i,,.e*'=“ + ^ 

Now, setting 


fc >0 


fc>3 


£1 =^£i,fce*'=“ + ^£i,_ 


fee 


—ike 


(134) 


k>0 


k>3 
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we can write 


3? 


{Lie*'^}=3?{( 


Li,_ie-*“ + Li _2e-2“+i;i 


)«“}■ 


By definitions (84) we have 

Ly_i = cf 7V(1, 2,1) + ^ cf-^’i7V(l, 1,1-1) 


l>3 


Ly_2 = cl’^N{l, 2,0) + ^ cf-'’'7V(l, 1 , 1 - 2 ) 


l>3 


If we set 


l>3 

E4 = J2cf~^'^N{l,l,l-2) 


(135) 

(136a) 

(136b) 

(137a) 

(137b) 


l>3 


we obtain just (133) from equations (136) and (1351. Once we have obtained the formula (1331, it 


only remains to bound properly the errors Ei, E^ and E 4 . From equation (134), by the triangle 
inequality and lemma we have 


< e 


fc >2 


l^'ll < 1^1,ol + 1^1,l| + 1-^1,fel + 1-^1,-fcl 

k>2 k>3 

,-g 3/3 Uoee^G-^/^ + 2^(1 + 

^ 6 ^ 2 Y^ 2 ‘^k^k-lQk-l /2 

3 

2^°ee^G-^^‘^ + 2^e(l + e)^G-'^^^ + 2®e(l + e)^G-^^^ + 2^^e^e^G^^^ 

eG-3/2 + G-"/2 ^ q-9/2 ^ g2g5/2 
gG-3/2 + q-7/2 ^ g2^5/2 


k>3 


< e-^ /3 

< 2i7e-G73 

< 2i®e-^'/3 


(138) 


Now we proceed with E 3 and E 4 . By propositions 12 and 16 from equations (137) 


|^^3| < ^|cf”^’^7V(l,Z,;- 1)1 < 23e^'^^ee"^"/^G^/^^(2‘^G"^)' < 2lVe■^"/^G■^/^ 

l>3 l>3 

\E 4 \ < ^|cf”^’^iV(l,l,;-2)| < 2ee'^^^ee-^"/3G^/2^(2^G-i)' < . 


l>3 


l>3 


The two estimates above, together with estimate (138) provide the desired bounds for the errors 
of equation (133). □ 


Putting together Lemmas and we already have 

£ = Lo + 23?{[(c?’^Af(l, 2,1) + £ 3 ) 6 "*“ + (c?’"iV(l, 2,0) + £ 4 ) 6 -"“ + £;i]e *®}+£>2 (139) 

with 

|£i (a, G; e) I < 2 IS (g'^/^ + e^G^/^ ^ ^^- 3 / 2 ^ ^- 0^3 

\E 3 {a,G-,e)\ < 220g"3/^ 

|^^4(Q;,G;e)| < 2^^eG^/^ 

\C> 2 {a,G,s-,e)\< 2 ^^G^/^e-^"l. 


(140) 
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We now compute N{1, 2,1) and N{1, 2,0). 


Lemma 28. Let N{q,m,n) be defined by equations (85(.Then 

Nil, 2 , 0) = + 2^^^ 


where 


\^Ett\ < 2®9G■2e■'^/^ \'^Ett\ < 2®9e"‘^/^ 

Proof. From proposition we have 

d2.1„_G3/3L2,l ^2.1 


Nil, 2 ,l) = ^e-^^/^ 


^4 WTT 


”^2^.1 “^2,1 


(141) 


where 

and 


< 45 2®G"3 |i?2il<18G3 


^(l>2,0) = '^e ^-22d2’°v^Y^+y'^o°v^(y^)^ + ^2,o + -Rlo 


(142) 


where 


ir^ipl < 45 26G-3 |i?i_o|<18G3. 

Taking the dominant terms in ( 141[ ) and (1421 we get: 

7V(1,2,1) = d2,ido’^^V^G-i/"e-^'/3 + ^Ep ^E 


(143) 


where 


and 


^E = 2-^d2,iV^{dl-^G-^ - 4’^G-5)e-^ ^Etr = (TI^ + Rij)d2,iG-^e-^ 
7V(1,2,0) = d2,odo’°^V^G3/2e-G'/3 + ^E + ^Etr 


(144) 


where 


^E = 2id2,ov^(d4’°G-S - 4’°G-i)e-^'/3 2 ^^^ ^ ^ R\ ^)d2,oG-^e-^"l\ 


Using the bounds given in Lemma 14 for dj ' and the bounds given in Lemma 13 for dm,n'- 

\^E\ < 2i|d2.i|v^(|di’"| + |4^|)G-ie-G'/3<27 9G-ie-«'/3 

\^Etr\ < M2.i|36G-2e-GV3<259G-2e-«'/3 


and also: 


\^E\ < 2i|d2,olv^(Mf| + M2’°l)G-3e-G^/3<269G-ie-GV3 

\‘^Etr\ < |d2,o|36e-G'/3 < 249e-G'/3. 


Using Lemma 14 d™’" = 1/(21)^"+^ and by definition (1061 for dm,n we have that 


^2,ic(q’ — —i 2 


d2.odo'^ — i‘2^ 


-l/2\ /-l/2\ / i 


2 

- 1/2 

2 


1 /V23/ 24 


2 ) 22 ' 
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We can then write equation (1431 as 


where 

satisfies 




^ Exx — ^ E ^ Ex 


\^Ett\ < 2^9G"5e"^"/^ + 2®9G"^e"^"/^ < 2®9G"^e"^"/^ 


In an analogous way, equation (144) can be written as 


iV(l, 2,0) = + "^Ett 


(145) 


(146) 


where 

satisfies 


2 77' _ 2 77’ I 2 77' 

-lL/TT — “T ^TR 


^^Ett] < 269G-5e-« /3 ^249e-^ < 2^'de-^ l^ 

and this proves the lemma. 

Using the approximations given in Lemma we have from lemmas |26| and \T7\ 


□ 


Lemma 29. For G > 32 and eG < 1/8, the Melnikov potential C (831 is given by 


£ = Lo + cos(s - a) + £3 + £ 5 ^ 

+ cos(s - 2 a) + £4 + ^e) + 23fi{F;ie*®} + £>2 


where £>2 and E^ with fc = 1,3,4 are given in equations (140) and 


UsI < 2^3 9 G"^e"'^"/^ \Ee\ < 2^^ 9 eiLe"^"/^ 


Proof. By lemma 28 we have that A^(l,2,1) and iV(l,2,0) are real and then coincide with their 
real part. Equation (139) gives the correct estimation of £. To complete the proof is enough to 
take 

£5 = c^' ■ ^Ext and £g = cf * ‘^Exx 


where ^Exx and '^Exx are given in lemma 28 Therefore by proposition 12 we find directly the 
bounds of £5 and Eq. 

□ 


Lemma 30. For G > 32 and eG < 1/8, the Melnikov potential £ (83) is given by 


£ = £, 


+ cos(s - a) + £3 + £5 + £7^ 

- cos(s - 2 a) (3v^eG3/2e-‘^"/3 + £4 + £5 + £ 3 ) + 23fJ{£ie*®} + £>2 


where £>2 and £fc with k = 1, 3,... 6 are given in equations (140) and 


|£7| < 98e2G-i/2e-GV3 < q^2^^^q3/2^-gV3 
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Proof. From lemma we have 

cf ^ y|g-i/2g-G3/3 ^ 


with 


Es = Q2V^G3/2e-GV3 


Therefore by lemma (291 and the bounds of Qi and Q 2 given in lemma 23 we conclude the proof. □ 


It only remains to estimate the Fourier coefficient Lq = £q defined in (55) or (831. 


Lemma 31. Let N{q,m,n) be defined by equations (85). Then for m,n S N, m + n > 0, 

\N{0,m,n)\ < 2 "^+"+ 2 g- 2 ™- 2 "+i^ 


Proof. Since r S K in the integral (85), it is easy to see that 

1 


1 


< 1 


r + z \T-i\ 


and then 


1 


1 


r + ip” r — i 


< 


1 


1 + r2 


For n,m > 0, using equation (85) and lemma 13 to bound dm,m the lemma follows: 

dr 


/ ° 

-C 


/ 1 H” 

/ —00 ~ ' 

_ 2n+lg—1/2^ ^ 2m—2n+l 


□ 


Lemma 32. Let k gN and Lo,fc defined by equation (55 K Then 

Lo,k= Y. Co'-"’-"A^(0,Z-fc,0 


l>k+l 


Proof. From equations (84), we have just to prove that for k > 2 

iV(0,0,fc) = iV(0,fc,0) = 0. 


By equations (85) this reduces to show that 

dr 


/_oo (t ± i) 


2 k 


= 0 


where the positive sign in the denominator correspond to /(O, 0, fc) and the negative to l{0,k,0). 
Since the variable t G K this integral is trivial 


dr 


1 


/_oo {t i i)2^ 2fc — 1 (t ± i)2fc-l 


= 0 


this proves the lemma. 


□ 
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Lemma 33. Let Lo.fe be defined by equations (841 for k>0. If G > 32, 

\Lo,k\ < 22'=+8e'=G-2'=-3. 

Proof. From lemma we have 

\LoM< E \co~^’~''\\N{0,l-k,l)\, 


and by propositions 12 and 16 


\Lo.±k\ < 2-2fc+3e'=G2'=+i 


-2k-3 


l>k+l 


□ 


Lemma 34. Let Lq = Cq be defined by equations (55) or (|83[). Then for G > 32 


Lq — Lo,o + tT^'G ^ + F 2 ) cos(a) + F 3 


^0,0 — Cq —Lt 


Fi 


where 


\Fi\<2^‘^G-\ IF 2 I < 2i*5eG-®, 

Proof. From proposition we know that 

Lq = Lq.o + 2 E -bo.fc cos ka, 


F3I < 2iVG-^. 


k>l 


and from lemma [32] we have that 


,0,0 = cyiv(o, 1,1)+ E co'’° w 0 0 

l>2 

(147a) 

, 0,1 = cg’-'iV(0,1, 2) + E cE'’”'^(0, ^ - 1, 0 

l>3 

(147b) 

, 0 fe= E Cq~’"'~'"N{0 ,1- kj) fork > 2 . 

(147c) 


l>k-\-l 


Introducing 

Fi = Eco'’°^( 0 ,^ 0 , Fa = 2 ^cE^’”^^( 0 ,^- 1 , 0 , F 3 = 2 ^ cosfcaLo,fc, 


l>2 


l>3 


k>2 


by lemmas 31 33 and proposition |12[ using the hypothesis on G we have 

|Fi| < 23 gE(2'‘G-‘‘)' < 212^--^ 


l>2 


IF 2 I < 2^eG^'^{2'^G-y < 2^^eG 

l>3 

IF3I <2^|Lo,fc| <2iVG-0 


k>2 


Now, from definition (851 we have that 


iV(0,l,l)= — 


2^ /-l/2\ /-1/2 


dr 


iV(0,l,2) = 


G^v 1 j\ 1 ;7-oo(^" + i) 

23 /-l/2\ /-l/2\ f°° dr 


( b 

( b 

V 2/ 

V 2/ 


G5 V 1 




/_oo (t- 0 (t + 0 ^ 


From these equations, substituting equations (147) in the definition of Lq and the bounds given in 
equations (148) we have proven this lemma. □ 
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A refinement of this lemma is 


Lemma 35. Let Lq = Cq be defined by equations (55) or (83). Then if G > 2^^^ 

15 

Lo = Lo,o + + ^2 + A’a) cos(a) + F3 

O 

^0.0 = 2^ ^ + Fi + F4 
where Fi, F 2 and F 3 are given in lemma and 

IF 4 I < 298G■^e^ lAsI < 22 98G"®e^ 


Proof. In lemma 23 we have computed the constants and Cg’^, then by setting 

F4 = ^QzG~^, F 5 = ^7rQ4G“®cosa, 

and using the bounds for Qg and <54 we find the desired bound for F 4 and F 5 . 

With this lemma we can rewrite lemmaexactly as Theorem and so we have proved it. 


3.1 


□ 


Acknowledgments 

The authors are indebted to Marcel Guardia, Pau Martin, Regina Martinez and Carles Simo for 
helpful discussions. 


References 

[Arn64] V. I. Arnold. Instability of dynamical systems with several degrees of freedom. Soviet 
Mathematics, 5(5):581-585, 1964. 

[B 0 IO 6 ] S. Bolotin. Symbolic dynamics of almost collision orbits and skew products of symplectic 
maps. Nonlinearity, 19(9):2041-2063, 2006. 

[CZll] M. J. Capihski and P. Zgliczyhski. Transition tori in the planar restricted elliptic 
three-body problem. Nonlinearity, 24(5):1395-1432, 2011. 

[DGOO] A. Delshams and P. Gutierrez. Splitting potential and the Poincare-Melnikov method 
for whiskered tori in Hamiltonian systems. Journal of Nonlinear Science, 10:433-476, 
2000 . 

[DGR13] A. Delshams, M. Gidea and P. Roldan. Transition map and shadowing lemma for 
normally hyperbolic invariant manifolds. Discrete Contin. Dyn. Syst, 33(3):1089-1112, 
2013. 

[DLS06] A. Delshams, R. de la Llave and T. M. Seara. A geometric mechanism for diffusion 
in Hamiltonian systems overcoming the large gap problem: heuristics and rigorous 
verification on a model. Mem. Amer. Math. Soc., 179(844):viii-|-141, 2006. 

[DLS08] —. Geometric properties of the scattering map of a normally hyperbolic invariant 

manifold. Adv. Math., 217(3):1096-1153, 2008. 

[Erd56] A. Erdelyi. Asymptotic expansions. Dover Publications, Inc., New York, 1956. 

[FGKR14] J. Fejoz, M. Guardia, V. Kaloshin and P. Roldan. Kirkwood gaps and diffusion along 
mean motion resonances in the restricted planar three-body problem. To Appear in J. 
Eur. Math. Soc., 2014. 


40 





[GL06] 

[GLS14] 

[GMS12] 

[GMS15] 

[LMS85] 

[LS80] 

[McG73] 

[Moe02] 

[Moe07] 

[MP94] 

[MS14] 

[Rob84] 

[Rob88] 

[Win41] 

[Xia92] 

[Xia93] 


M. Gidea and R. de la Llave. Topological methods in the instability problem of Hamil¬ 
tonian systems. Discrete Contin. Dyn. Syst., 14(2):295-328, 2006. 

M. Gidea, R. d. 1. Llave and T. M. Seara. A general mechanism of diffusion in hamil- 
tonian systems: Qualitative results, preprint, arXiv:1405.0866, 2014. 

M. Guardia, P. Marthn and T. M. Seara. Oscillatory motions for the restricted planar 
circular three-body problem, preprint, arXiv:1207.6531, 2012. 

—. A lambda lemma and oscillatory orbits for the restricted elliptic planar three 
body problem, preprint, 2015. 

J. Llibre, R. Martinez and C. Simo. Tranversality of the invariant manifolds associated 
to the Lyapunov family of periodic orbits near L 2 in the restricted three-body problem. 
J. Differential Equations, 58(1):104-156, 1985. 

J. Llibre and G. Simo. Oscillatory solutions in the planar restricted three-body problem. 
Math. Ann., 248(2):153-184, 1980. 

R. McGehee. A stable manifold theorem for degenerate fixed points with applications 
to celestial mechanics. J. Differential Equations, 14:70-88, 1973. 

R. Moeckel. Generic drift on Gantor sets of annuli. In Celestial mechanics (Evanston, 
IL, 1999), volume 292 of Contemp. Math., pages 163-171. Amer. Math. Soc., Provi¬ 
dence, RI, 2002. 

—. Symbolic dynamics in the planar three-body problem. Regul. Chaotic Dyn., 
12(5):449-475, 2007. 

R. Martinez and C. Pinyol. Parabolic orbits in the elliptic restricted three body problem. 
J. Differential Equations, lll(2):299-339, 1994. 

R. Martinez and G. Simo. Invariant manifolds at infinity of the RTBP and the bound¬ 
aries of bounded motion. Regul. Chaotic Dyn., 19(6):745-765, 2014. 

G. Robinson. Homoclinic orbits and oscillation for the planar three-body problem. J. 
Differential Equations, 52(3):356-377, 1984. 

—. Horseshoes for autonomous Hamiltonian systems using the Mel'nikov integral. 
Ergodic Theory Dynam. Systems, 8* (Charles Conley Memorial Issue):395-409, 1988. 

A. Wintner. The Analytical Eoundations of Celestial Mechanics. Princeton Mathemat¬ 
ical Series, v. 5. Princeton University Press, Princeton, N. J., 1941. 

Z. Xia. Melnikov method and transversal homoclinic points in the restricted three-body 
problem. J. Differential Equations, 96(1):170-184, 1992. 

—. Arnol'd diffusion in the elliptic restricted three-body problem. J. Dynam. Differ¬ 
ential Equations, 5(2):219-240, 1993. 


41 


